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Abstract 

, A sixth-order quadrupole boson Hamiltonian is used to describe the states 0"*" and 2+ identified 

I in several nuclei by various types of experiments. Two alternative descriptions of energy levels 

are proposed. One corresponds to a semi-classical approach of the model Hamiltonian while the 



other one provides the exact eigenvalues. Both procedures yield close formulas for energies. The 
first procedure involves four parameters, while the second involves a compact formula with five 
parameters. In each case the parameters are fixed by a least-square fit procedure. Applications are 



00 
O 

> 

■ performed for eight even-even nuclei. Both methods yield results which are in a surprisingly good 

agreement with the experimental data. We give also our predicted reduced transition probabilities 
within the two approaches, although the corresponding experimental data are not yet available. 

PACS numbers: : 21.10.Re, 23.20.Lv, 21.60. Ev 
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I. INTRODUCTION 



The collective states of deformed nuclei are usually classified in rotational bands dis- 
tinguished by a quantum number K, which is the angular momentum projection on the z 
axis of the intrinsic reference frame. The collective character of the states is diminished by 
increasing the value of K [l|,|2|,y,3. In Ref. [si one of us (A.A.R.) suggested a possible 
method of developing bands in a horizontal fashion. Indeed, therein on the top of each 
state in the ground band a full band of monopole multi-phonon states has been constructed. 
The states of the newly constructed band on the top of the ground band state of angular 
momentum J have the same angular momentum J. This feature contrasts the property 
of the ground band, where the states have different angular momenta, i.e. 0, 2, 4, 6 etc. 
This idea has been recently considered in a phenomenological context trying to organize 
the states, describing the motion of the intrinsic degrees of freedom, in bands. Thus, two 
intrinsic collective coordinates, similar to the nuclear deformations /3 and 7, are described 
by the irreducible representations of a SU(2) group acting in a fictitious space (i.e. not in 



ordinary space) . Compact formulas for the excitation energies have been obtained 

Recently, about 26 states 0^ and 67 states 2+ have been populated in ^^^Er by means of a 
{p,t) reaction [Sj. In the cited paper the excitation energies and the corresponding reaction 
strength have been provided. These data were described qualitatively by two microscopic 
models, called projected shell model (PSM) and quasiparticle phonon model (QPM), respec- 
tively. Both models have some inherent drawbacks. PSM restricts the fermion space to four 
quasiparticle states and even from the four qp space the states with four alike quasiparticles 
are excluded. This is not the case of QPM where the multi-quasiparticle components are 
taken into account by means of the QRPA approach. However, the final states contain at 
most two phonon states. These states violate the Pauli principle and moreover are not states 
of good angular momentum. 

In Ref. [9] some of us made a first attempt to fit the data of Ref. [8] using a phenomeno- 
logical model, namely a sixth-order quadrupole boson Hamiltonian that was developed in 
Ref. [7]. Since then about 12 new 2+ states have been identified by a more careful analysis 
of the data produced in the {p,t) experiment [10]. Here we show that the complete O"*" and 
2+ data sets, presently available, are nicely described by the closed formulas provided by 
the model of Ref. [7]. 
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Here we present details about both the semi-classical approach and the boson descrip- 
tion of these states. The compact analytical formulas are used to explain the data about 
the states 0+ and 2+ in several even-even nuclei: ^^^^^^^Gd, ^^^Dy, ^'^^Er, i^'^Hf, ^^^'^^^W, 
^^'^Os. The model quadrupole boson Hamiltonian is presented in Section 2. Therein we 
also present two distinct approaches for its spectrum. Analytical formulas for the reduced 
transition probabilities, corresponding to the mentioned treatments, are derived in Section 
III. Numerical applications to eight nuclei are presented in Section IV. The final conclusions 
are summarized in Section V. 



II. THE MODEL HAMILTONIAN 

We attempt to describe the set of states 0'*' and 2^ identified in various experiments in 
terms of quadrupole bosons, by means of the model Hamiltonian: 

H = eN+ Y: Cj[{blbl) {b,b2)j] + F {blbl) N {b,b2),, (2.1) 

J=0,2,4 

where b^^j,, &2^t! with —2 < /i < 2, are the quadrupole boson operators and N the boson 
number operator. The first remark about the chosen Hamiltonian refers to the fact that it 
commutes with the boson number operator. We recall that this feature is one of the signatures 
of the interacting boson approximation (IB A) [l^ which, as a matter of fact, was very 
successful in describing rotational bands in non-spherical nuclei. Moreover, the Hamiltonian 
given by Eq. (12.11) with F = 0, has been used, even before the IBA was proposed, to describe 



the yrast bands in transitional and deformed nuclei [20|, |2l| . Thus, an analytical formula for 
the yrast ener gies has been obtained, which in fact was generalizing the empirical expression 
used by Ejiri 22]. As in Ref.Q], this Hamiltonian is alternatively treated semi-classically 



and exactly solved in the boson space. For a self-contained presentation we give here the 
basic results obtained in the mentioned treatments. 

A. Semi-classical treatment 

The boson Hamiltonian (12.11) is treated by a Time Dependent Variational Principle 
(TDVP): 

6j\^\ (^H-th-^^ \^)dt' = 0. (2.2) 



When the variational state spans the whole space of the boson states, solving the equa- 
tion (12.21) is equivalent to solving the time dependent equation associated to the model 
Hamiltonian H. The classical features encountered by H can be described by restricting the 
space of 1^) to the coherent states: 

1^) = exp [zobS - z*bo + Z2(b5 + hU) - z;(b2 + b_2)] |0). (2.3) 

Here the boson vacuum state is denoted by |0). The function depends on the complex 
parameters zq,Z2 and their complex conjugates Zq,Z2. These parameters play the role of 
classical phase space coordinates whose equations of motion are provided by the TDVP 
equations. By a suitable change of coordinates, 



= 2('=+2)/^/?e(zfc), p. = h2^^+^^'Hm{z,), 
k + 2 

k = 9, i = — — , 



(2.4) 



the classical equations of motion acquire a canonical form, while the classical Hamilton 
function (the average of H with I^P)), 7i, becomes a function of the generalized phase space 
coordinates, q and p: 



1 



C 



+ 
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4 + (l2 + ^ipl+ pDJ - ^(^1^2 - q2Pl] 



ql + ql + j^ivl+vl) 



(2.5) 



where the factors A, B and C have simple expressions in terms of the coefficients e, Cj 
involved in the boson Hamiltonian: 



A = e, B = -Cq + -^C2 + — C4, 
5 " 7^5 35 ^' 
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^ - ~5^° + ~ 35^^- 



(2.6) 



Conventionally, we shall call the part of Ti not depending on momenta, as the potential 
energy of the system: 



V{q) = n\p,=p^=o, with (q) = (qi,q2). 



(2.7) 
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Thus, the potential energy associated to Ti is: 

V{q) = I (g? + ql) + f (g? + qlf + ^ (g? + • (2.8) 

In Ref. j3] we showed that this classical function exhibits a symmetry with respect to the 
classical rotations generated by the classical functions obtained by averaging the generators 
of a SUbi2) algebra with |*): 



L3 



n r 
4 - 

h 



2V2 
n 



2^/2'. 



2h%-{h\ + hU){h2 + h^2 
bl{b2 + b_2) + ibl + bU)bo 

'blib2 + b^2)-ibl + bU)b, 



(2.9) 



Thus, the generators of the classical SUc{2) algebra acting in a fictitious space are defined 
by: 

A. = (^|Lfc|^),A; = l,2,3. (2.10) 

It can be checked that the classical system has two constants of motion and these are 7i and 
£3. On the other hand the system is fully described by two degrees of freedom, qi and q2. 
Consequently, the classical system is fully solvable (or integrable). Therefore, the equations 
of motion can be integrated and the trajectories analytically described. 

7i contains two distinct terms describing an anharmonic motion of a classical plane oscil- 
lator and a pseudo-rotation around an axis perpendicular to the oscillator plane, respectively. 
Taking into account that the third component of the pseudo-angular momentum is a con- 
stant of motion, the classical Hamiltonian considered in the reduced space can be easily 
quantized and the resulting energy is: 



en,M = A{n + 1) + B{n + if + + ^ 

2 5 



[n + l] 



4(n + 1)M^ 



(2.11) 



The number of the oscillator quanta in the gi, g2 plane is denoted by n while the value of 
the third component of the pseudo-angular momentum is M. Actually, Eq. (12.111) represents 
a semi-classical spectrum which describes the motion of the intrinsic degrees of freedom qi 
and g2, related to the nuclear deformations [3 and 7. 

Assuming that the rotational degrees of freedom are only weakly coupled to the motion of 
the intrinsic coordinates, the total energy associated to the motion in the laboratory frame 



can be written as a sum of two terms corresponding to the intrinsic and rotational motions, 
respectively: 

en,Af,j = A{n + 1) + B{n + 1)^ + ^M^ 



+ 5 



n+lf -4:{n + l)M'^ +SJ{J+1) (2.12) 



Averaging both angular momenta squared, and L^, on one obtains a relationship 



between the two momenta. Thus (see Ref. 



7|), to the values J = and J = 2 correspond 



different values of M, namely M = and M = 1, respectively. 

Therefore, considering the above equation for the sets of states with angular momenta 
J=0, 2 and normalizing the results to the energy of the first 0"'' state, one obtains the 
following expressions for the excitation energies: 

Enfi = ^n,o,o - ^0,0.0 = Ifu' + (|f + + {A + 2B + |F)n, n > 0, 

5 5 5 

En,i = K,i,2 - ^0,0,0 = Ifu' + (|f + By + (A + 2B- ^F)/! + C, n>l, (2.13) 

5 5 5 

where 

B. Exact eigenvalues 

Note that the model Hamiltonian is highly anharmonic due to the terms of fourth and 
sixth-order in the quadrupole phenomenological bosons. Despite this fact it is easy to see 
that this Hamiltonian is diagonal in the boson basis \NvaJM) , where the quantum numbers 
have the significance of the boson number (iV), seniority (f ), missing quantum number (a), 
angular momentum (J) and its projection on the axis OZ (M). These basis states have been 
analytically studied in Ref. jlj] using alternatively different representations like, laboratory 
frame coordinates, intrinsic frame coordinates, boson variables. To prove the statement 
concerning the diagonal form of H in the mentioned boson basis, it is useful to write the 



fourth order term in a different form (see Ref. 
expression for H: 



20l |) which results in having a more convenient 



H = (A + 7)iV + (5 + ^)iV2-i(i? + ^+7) J' 

- (44)jM2)o + i^ (44), (6262)0, (2.15) 



where the coefficient 7 has the expression: 

7 = ^^2-1^4. (2.16) 
^ 7^5 T ^ ^ ^ 

From this expression it is obvious the H commutes with the operators N, A, J^, where 
A denotes the Casimir operator of the group R^: 

A = N{N + 3)- 5{blblUhh)o. (2.17) 

The eigenvalue corresponding to the state \NvaJM) is: 

En,v,j = \fN' + (5 + ^F)N^ (2.18) 

+ (A + 7 - k\c + \f))n - 1(5 + + 7) J( J + 1) 

+ {\c + + 3(ic + \f)v - \fNv^ - \fNv. 

8 5 8 5 5 5 

Comparing this with Eq. fl2.12p . we notice that the eigenvalues of if, corresponding to a given 
J, are characterized by two quantum numbers, namely the number of bosons A^ and the 
seniority v. Therefore, using the new expression for energies one expects a better description 
of the data. For J = we use the lowest two values for seniority quantum number, i.e. 
= 0, 3, and obtain: 

ENfifl = \fN' + {B+ ^-F)N^ + (A + 7 - - ^-F)N, AT = 0, 2, 4, ... (2.19) 
i^iv,3,o = \fN' + {B + ^F)N^ + (A + 7 - - '^F)N + + ^F, AT = 3, 5, 7, .... 

Similarly, for J = 2 we consider the lowest two allowed seniorities, i.e. f = 1, 2. The result 
is: 

En,i,2 = ^FAr3 + (5 + ^F)Ar2 + (A + 7-^C-2F)Ar-5-7 + ^C + ^F,Ar = l,3,5,... 

11 3 16 9 

Enx2 = -FN' + {B + -F)N^ + {A + - -C - —F)N - B - + -C + AF,N = 2A.Q. ■■■■ 
5 5 8 5 8 

(2.20) 

Note that within the boson treatment the lowest energy denoted by -Eo,o,o is equal to zero 
and therefore there is no need to renormalize the energies of excited states with respect to 
the ground state energy. Let us now turn our attention to the missing quantum number a. 
This quantum number labels the irreducible representations (J) which appear in an R^ 



irreducible representation (v). The name is suggesting that there is no intermediate group 
between and R3 whose Casimir operator might make the distinction between different J 
representations corresponding to the same seniority v. a labels the solutions of the double 
inequality for the integer number p : 



The number of solutions for this double inequality is the degeneracy d^{I), characterizing 
the reduction D R3. It is clear that for J = 2 and fixed v, the number of solutions of Eq. 
(2.21) is either or 1. For example, there is no state 2"*" with v = 0,3,6,9, .... Concerning 
the states J = 0, one has dy{0) = 1 if v = 3k with k positive integer and dv{0) = otherwise. 

Concluding, for J = 0, 2 there is no degeneracy, i.e. the set (f , J) either does not exist 
or is uniquely determined by the relation fl2.2ip . 

III. ELECTRIC QUADRUPOLE TRANSITIONS 

The states 0'*' can be related to the states 2^ by E2 transitions whereas the states of 
the same angular momentum are related by EO transitions. Since the EO transitions for 
highly excited states are not yet experimentally investigated we confine our study to the 
E2 transitions. Compact formulas for E2 transitions have been presented in our previous 
publication [9]. However details about the derivation of these expressions were not given. 
Here we complete the description of the E2 transitions by providing additional information 
which will facilitate a straightforward derivation of the results listed in the reference quoted 
above. 

A. Semi-classical approach 

We suppose that the leading contribution to the E2 transitions is provided by the linear 
boson term: 



The average of this operator with the coherent state |\E') (see Eq. (l2.3l) )has the expression: 



V — J < 3p < V 



, for J = even, p = integer. 



(2.21) 




(3.1) 



(3.2) 
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Here 5m,n stands for the Kronecker symboL 

The semi-classical energies have been obtained by quantizing the plane oscillator defined 
with the coordinates qi and q2- Thus, the energies depend on the total number of quanta 
along the two plane axes. It is convenient to use the polar coordinates associated to the 
Cartesian qi and q2- The principal and radial quantum numbers are related by: 

2nr + (5j,2 = n. (3.3) 

Since the M-quantum number is equal to for O"*" states and 1 for the states 2"*", one can 
use only one label for the intrinsic states jn^, M): 

|0n) = |^,0), |2.) = |^,1). (3.4) 

Using the explicit wave functions for the plane oscillator one calculates the matrix elements 
of the function 7^^. In the laboratory frame, the transition operator is acting on both the 
coordinates gi, q2 and the Euler angles VL = {Pi,62,9^) and has the expression: 

T2M = g.v^ [qiDU + ^{Dln + Dl^_^)^ , (3.5) 

where Dj^jj^ is the Wigner function describing the rotation matrix. 

In the liquid drop model the state of angular momentum 2 in the laboratory frame 
consists of two factors, one depending only on the deformation /3, while the other one is 
linear combination of the Wigner functions -Df/^ with the coefficients qk depending on 
the deformation 7. In the present formalism by averaging the boson Hamiltonian on the 
coherent state one obtains the equations of motion for the intrinsic variable gi, q2 which 
may be related to the deformations (3, 7. Therefore, we assume that in the laboratory frame 
the wave functions are factorized in the following manner: 



\2nM) = |2„)^2A/(1^), ^2Mm = \h:^ 

V O ZTT 

lOriO) = |0„)^. (3.6) 



D\m + D\j2 + D\t_2 



Using the convention of Rose [17| for the reduced matrix elements, we have: 

B{E2-X' ^ 0+) = (2n'||T2||0n)2. (3.7) 

with n>l. Analytical expressions for the above B(E2) values as well as for some particular 
branching ratios were given in Ref. ^. 
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For a transition operator having an harmonic structure, the E2 transition between any 
two states 2"*" is forbidden. This result is specific to the present semi-classical description. 
Indeed, the matrix elements of the variables qi and q2 between the states |2„) and |2„') are 
equal to zero due to the integration over the polar angle. In Ref. tI we gave a group theory 
argument for this result. Indeed, with respect to the pseudo- rotation group the harmonic 
transition operator is a tensor of rank 1/2 while the states 2"*" have the pseudo-angular 
momentum equal to 1. Then, it becomes manifest that two states 2+ cannot be linked by 
an harmonic transition operator. Of course, that is not true in the boson treatment, as we 
shall see in the next subsection. 

In order to get non-vanishing transition matrix elements between two different 2+ states 
we introduced an anharmonic term in the expression of the transition operator: 



rpanh 

^2^l ~ Hanh 



(3.8) 



The reduced matrix element between two 2+ states can be analytically obtained [9| A 
peculiar feature of the present formalism is the fact that the anharmonic term does not 
contribute to the transition 2+ 0+. On the other hand, as we have already mentioned, 
the harmonic term does not contribute to the transition 2+ 2^,. Thus, the final result for 
the transition 2,^ 2^_2 is proportional to {n — Vf. 



B. E2 transitions within the boson picture 



In what follows we shall identify the missing quantum number with the integer positive 



number p which satisfies the inequality, 
\n\pIM) have a factorized form isl . 



fl2.2ip . In the intrinsic frame of reference, the states 



\nXpIM) = Fnx{P)gxpi{7,n), 



(3.9) 



where 



2(2(^-A))! 



f3"L'^^L)il3')expi-b'). 



(3.10) 



L^~^^_ , (/3^) stands for the generalized Laguerre polynomial. The functions Fn\ are orthonor- 



L(n-A) 



malized on the interval [0, oo) with the integration measure f3^df3. We need the normalized 
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unctions depending on the variables 7 and Q for the angular momenta and 2. These are 



13|, 



^000(7,^) 
^102 (7>^) 



1 


^310(7,1^) = 


471^2' 


1 /5 




47rV2 


cos7L'Afo + ■ 


1 /5 




47rV2 


cos 27L)^^o - 



1_ 3 
47rV 2 

sin 7 (j 
V2 ^ 
sin 27 



cos 37, 



)2 



{Din + D 



M,-2 



In the intrinsic frame, the harmonic transition operator has the expression: 



(3.11) 



qh(3 cos 71)^/0 + 



2 , sin7.n2 



(3.12) 



The reduced matrix elements of the transition operator between the states described in the 
previous section are calculated in Appendix A. The reduced probability for the transition 
\nvpIM) \n'v'p'I' M') is obtained by squaringthe corresponding reduced matrix element 



of the transition operator. As shown in Ref. 
simple. 



3] the final analytical expressions are very 



IV. NUMERICAL RESULTS 

The equations derived in the semi-classical framework fl2.13p for the energies of the states 
0"*" and 2+ were used to fit by a least square procedure the data for several nuclei: ^^^'^^^Gd, 
i62j3y^ les^r, 176 Hf, i80,i84^^ igoQs. The parameters, A, B, C, F yielded by the fitting proce- 
dure are listed in Table I. The fitting procedure provides also the set of quantum numbers 
{rifcjfc associated to the states specified by the ordering index k. Of course the values of n 
provided by the equations expressing the condition that the value is minimum are not 
integers. We assigned to a given energy level k the integer which is closest to yielded by 
the least square equations. 

The boson description provides the expressions fl2.19p for energies of the states 0^, while 
for the states 2~^, Eq. (12.201) is determining the energies. These equations define four sets of 
energies which are depending on five parameters: A, B,C, F,'~f. These parameters together 
with the quantum number N are to be fixed by a least square procedure. For comparison we 
performed the fitting procedure for the same nuclei considered in the semi-classical approach. 
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Amazingly, both procedures lead to a cubic expression in n and respectively, although the 
two quantum numbers have different significance. Indeed, the quantum number n represents 
the number of the plane oscillator quanta associated to the intrinsic degrees of freedom gi, q2 
related to the nuclear deformations (3, 7. On the other hand the quantum number is the 
number of the quadrupole bosons which are describing the system in the laboratory frame. 

Note that in principle both the boson number and seniority could be obtained by solving 
the least square equations but the procedure would be quite tedious. For the sake of sim- 
plicity we kept only the boson number as variable to be determined and chose the lowest 
seniorities. The reason is that for these seniority values the energy equation has a similar 
structure as in the semiclassical case. For the states 0"'' we started with the f = expressions 
and tried to describe all energies as corresponding to f = 0. The result was that a set of 
calculated energies exhibit large deviation from the experimental data. These states were 
considered to have v=3. In the next step both expressions, corresponding to v=0 and v=3, 
have been used with the assignments determined before, and new least square equations 
have been written for the five parameters and the A^s( number of states) values of A^ (boson 
number). 

The results of the fitting procedures concerning the structure coefficients mentioned above 
are collected in Table I . Inserting the fixed coefficients in the equations defining the energies, 
one obtains two sets of energies for classical description and four sets for the exact treatment. 
The sets of energies are plotted as functions of n and A^ respectively in Figs. 1-16. 

The four sets of energies, E^^^o, En^^^, -Eat,!, 2, En,2,2 with the restrictions for A^ mentioned 
above (see Eqs. 12.191 and I2.20p . are represented in the panels b), c) of the left figures and 
panels b), c) of the right figures, respectively. In the panels a) of the left and right figures 
the semi-classical energies for the states 0"^ and 2+ are given. The full line curve is the 
energy as function of A^, with A^ considered as a continuous variable. The integer number 
which lies closest to the experimental data is the assigned quantum number A^. We remark 
that the agreement with the experimental data is quite good for both semi-classical and the 
exact eigenvalues. The remarkable feature of our approach is that by compact formulas we 
obtain a realistic description of a large number of excitation energies, despite the fact that 
the number of the fitting parameters is relatively small. 

The states studied in this paper have been populated, by several groups, in experiments 
like (p, t), it,p), {D,D'), (n, 7), {n,n'), Coulomb excitation. 
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-21.281 
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343.827 


-20.827 
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294.746 
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-40. 
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31.978 
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621.2 
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-280. 
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83.2319 
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24. 


0.3865 


299.8638 
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-111.787 
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841.961 


180^ 


1030.37 
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-656. 


11.45 


104.28 


-87.74 


-40 


11.45 


747.06 




781.966 


-50.0314 


-523. 


4.3282 


111.469 


-48.3001 
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4.3282 


596.225 


i90Os 


856.71 


-52.904 


-567.372 


4.51 


133.332 


-50.472 


-324. 


4.51 


504.188 



TABLE I: The structure coefficients yielded by the fitting procedure apphed to the semi-classical 
expressions for the 0^ and 2+ energies are given in the first four columns while the results for those 
involved in the exact eigenvalue expressions are given in the last five columns. All coefficients are 
given in units of keV. 



For ^^^Gd, the energies of 11 states O"*" and 34 states 2"*" are known from Ref. 2J, |25|. The 
energies for 0"*" states are smaller than SOOOkeV while the states 2+ lie below 3500keV. In 



24 



27 



28 



29 



30|. The results 



the case of ^^^Gd, 15 O"*" and 60 2+ are known from Refs. 
for ^^^Gd are plotted in Figs. 1,2 while those for ^^^Gd are represented in Figs. 3, 4. 

Classical results can be interpreted in terms of quantized states of classical trajecto- 
ries describing the motion in the potential V(q), defined by Eq. fl2.8p . This potential has 
been plotted in Fig. 17 for ^^^Gd and ^^"^Gd, respectively and in Fig. 18 for ^^^Er, using 
alternatively the set of parameters provided by the semiclassical description and the exact 
treatment. 

We recall the fact that along the isotopic chain of Gd one records a transition from 
spherical nuclei (the light ones) to deformed like nuclei. The first set of nuclei satisfy an 
SU{5) symmetry while the second one an SU{3) symmetry. The critical nucleus for this 
transition is considered to be ^^^Gd, which itself exhibits a distinct symmetry called X{5) 
symmetry. This transition critical point is characterized by a specific value for the ratio 
E4+ / E2+ and special features in the E2 properties of the ground as well as of the adjacent 
bands. The question is whether we find some fingerprints for this shape transition in the 
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FIG. 1: Excitation energies of the J'^ 



FIG. 2: Excitation energies of the 



0+ states, in ^^^Gd, described semi-classically, states, described semi-classically, panel a), and 

panel a) , and by eigenvalues of the model Hamil- by eigenvalues of the model Hamiltonian, cor- 

tonian, corresponding to the seniority v=0, responding to the seniority v=l, panel b) and 

panel b), and v=3, panel c), are compared with v=2, panel c), are compared with the experi- 



the experimental data. 



mental data. 



semiclassical description. This is in fact the reason we present here the potential energy 
corresponding to the two even isotopes of Gd. 

Comparing the potential for the two Gd isotopes considered here, we note that the slope 
of V{q) in its ascending part is higher for ^^^Gd than for ^^^Gd. In other words the first 
minimum is more flat for ^^^Gd than for ^^^Gd. The consequence is that the first 2"^ state 
is lower in energy for ^^"^Gd than for ^^^Gd. Since the transition for the state 4+ is felt less 
strongly than in the state 2"*" 23|], the ratio of the two states energies is of course seriously 
affected. The second remark refers to the fact that the secondary minimum for ^^^Gd is 
more deformed than that corresponding to ^^^Gd. 

In ^^^Dy have been observed 12 states O"*" and 11 states 2" 
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3l[ |. The results of 



our calculations are compared with the corresponding experimental data in Figs. 5, 6. We 
note that by contrast to the case of Gd' s isotopes here the energies are distributed around 
the maximum of the curves. Also we remark the low states density around the secondary 



14 



4000 




4000 




20 40 n 60 80 100 
^ 4000n 




20 40 M 60 80 100 



FIG. 3: The same as in Fig.l but for ^^*Gd. FIG. 4: The same as in Fig. 2 but for ^^*Gd 
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FIG. 5: The same as in Fig. 1 but for ^"^"'Dy. FIG. 6: The same as in Fig. 2 but for ^'"'Dy. 



minimum. 

In ^^^Er, by a {p,t) experiment, there have been identified 26 energy levels O"*" and 79 
states 2+ 

As in any other theoretical model, the number of predicted energy levels in our model is 
equal to the number of the considered basis states which, in general, is different from the 
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FIG. 7: The same as in Fig. 1 but for ^''^Er. FIG. 8: The same as in Fig. 2 but for ^''^Er. 

number of the experimentally identified states. Clearly, in the plots shown here there are 
more predicted energy levels than experimental energies. 

The predicted energies which do not have correspondent data may require higher resolu- 



tion or a different type of experiment. For example for ^^^Gd, the authors of Ref. 10| found 



several new states through an {p,t) experiment, that could not be seen by the previous 



(n, n') experiment 



Concerning the predictive power of the present formalism it is worth mentioning an 
interesting story concerning the case of ^^^Er. Indeed, after the publication of data in 
Ref. Q], where only 67 levels 2"*" have been reported, and shortly after we provided a 
phenomenological interpretation in Ref.^], Bucurescu and his collaborators analyzed more 
carefully the data and found another 12 energy levels with angular momentum 2 and positive 
parity 12|. Of course, this was a challenge for us since explaining the new data is indeed a 
severe test for the proposed theoretical description. These 12 new levels are also considered 
here, keeping the fitted parameters from Ref. jsl unchanged. 

The new data for the 2"*" energies are given on the first column of Table II. The first 
three values fill the vacancies in the curves of Ref. jol and are presented here in Fig. 8. 
The assigned quantum numbers A^, v are those given in the second column of Table II. The 
remaining data lie very closely to the data which are already represented in Fig. 8. The later 
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New energies [keV] 


N 


V 


Old energies [keV] 


(2174.0) 


63 


1 




2580.4 


70 


2 




2683.2 


71 


1 




2969.3 


47 


1 


2.961.2 


3391.1 


39 


1 


3361.9 


3418.2 


17 


1 


3429.2 


3794.1 


31 


1 


3789.5 


3838.0 


29 


1 


3861.9 


(3923.4) 


25 


1 


3933 


(4009.6) 


26 


2 


4033.5 


4060.7 


80 


2 


4055 


4069.2 


81 


1 


4075.6 



TABLE II: Twelve data for the newly identified 2^*" states are ^iven on the left column. The first 
three are falling on the graph representing the data from Ref.[S] and are interpreted as v=l and 
v=2 states of large N. The remaining energies from the first column lie closely to the old energies 



from Ref. 



described as low seniority states. 



are given also in the fourth column of Table II, together with the correspondingly assigned 
and v values. The quasi- degeneracy for the energy levels between 2.9 and 4.1 MeV, shown 
in Table II, may suggest t 
set of states \NvaJM > 
led us to the conclusion that the new energies from the first column might correspond to 
(A^, v) values which are different than those given on columns 2 and 3. This suspicion is 



lat a symmetry exists. However, as mentioned already before, the 
13, Q does not comprise any degeneracy for J = 2. This feature 



based on the nonlinear character of the equation in N and v 

En,vJ = 



(4.1) 



for a given value of £. Keeping the same parameters as before, we obtained the theoretical 
values for energies given in Table III. Concluding, the compact formula given by Eq.(2.20) 
may describe a large amount of data despite the fact that only few parameters are involved. 
Most of the data are described as low seniority states but for ^^^Er, there are also energy 



17 



New energies [keV] 


Theory 


N 


V 


2969.3 


2961.4 


45 


5 


3391.1 


3408.4 


39 


7 


3418.2 


3426.8 


14 


10 


3794.1 


3783.7 


17 


11 


3838.0 


3833.1 


80 


8 


(3923.4) 


3921.9 


28 


4 


(4009.6) 


4007.5 


26 


8 


4060.7 


4055.6 


26 


10 


4069.2 


4068.8 


27 


11 



TABLE III: The excitation energies for the newly identified 2"*" states, first column, are compared 
with the predictions of the energy expression (2.19), (2.20), given in second column, for N and v 
from third and fourth column, respectively. 

levels which correspond to high seniority states. 

Finally, we remark that the least square procedure yields for the first excited 0"*" state 
in ^^^Er, a value for the boson number equal to three. On the other hand in Ref.|l6], by 
means of a {n, 7) reaction a complete scheme of levels has been produced for J < 6 and 
< 2 MeV. The result is that there is no state 0^ with an energy smaller than 1.217 
MeV. Thus, it is an open question which deserves further consideration, whether there are 
specific selection rules which prevent the population of the predicted = 2, t> = state by 
the experiments mentioned above. 



In ^^^Hf , eleven 0"'' and seven 2+ states are available 2J, |32| . The theoretical results are 
plotted in Figs. 9,10 together with the experimental data. From there one notices that only 



few states of seniority 2 (1) and three 



For ^^°W one knows, from Ref. 
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2) are found. 



33| . the energies of eleven O"*" and five 2"*". Here we 
assigned the seniority 2 only to one state 2"*" and seniority 3 only to two states 0^. 

For ^^^W more states are experimentally known. Indeed, in Refs. 2J, |3J] the energies of 
sixteen O"*" and eighteen 2^*" have been reported. 

The last nucleus investigated is ^^''Os for which we know seven energy levels 0+ and 



seventeen 2"*". Experimental data are those from Ref. 
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34|. 
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FIG. 9: The same as in Fig. 1, but for ^^^Hf. FIG. 10: The same as in Fig. 2, but for ^'^^Hf. 




FIG. 11: Fig. 11. The same as in Fig. 1, but FIG. 12: The same as in Fig. 2, but for ^'^"W. 
for 180W. 
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FIG. 17: The classical potentials (2.7) corre- 
sponding to the sets of parameters specified in 
Table I for ^^^Gd (upper panel) and ^^^Gd (lower 
panel) provided by the semi-classical description 
are plotted as functions of the polar coordinate 



FIG. 18: The classical potentials (2.7) corre- 
sponding to the sets of parameters specified 
in Table I for ^^^Er, provided by the semi- 
classical (upper panel) and exact descriptions 
(lower panel) arc plotted as functions of the po- 
lar coordinate r(= \/{q\ + Q^))- 
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V. CONCLUSION 



In the previous sections, we proposed two phenomenological descriptions of the excitation 
energies of the states 0^ and 2"*" experimentally identified in several even-even nuclei. They 
correspond to two distinct ways of treating the same sixth-order quadrupole boson Hamil- 
tonian. One is a semi-classical description while the second one uses the exact eigenvalues. 
While in the yrast band the highest seniority states are the best candidates for a realistic 
description, for the states of the same angular momentum, the lowest seniority states are 
used for most states. We found, however, that some of the states 2^ of ^^^Er are higher 
seniority states. It is remarkable that both 0"'' and 2^^ states exhibit a cubic n dependence. 
We know that such a behavior for energy in the yrast bands is determining a back-bending 
{isl phenomenon for the moment of inertia as a function of the rotational frequency. Here 
a back-bending also shows up but the cause is different from that determining the bending 
in the moment of inertia in the yrast band. 

The terms of the classical Hamiltonian which do not depend on momenta define the 
potential of the classical system. This has been plotted in Fig. 17 for ^^^'^^^Gd and Fig. 18 for 
i68g],_ For Gd isotopes we used the parameters provided by the semi-classical treatment while 
for ^^^Er the two panels correspond to two sets of parameters obtained by classical and exact 
descriptions, respectively. From the upper panel of Fig. 18 we notice that some semi-classical 
states may accommodate the second well of the potential [7]. The boson description yields 
a similar spectrum as the semi-classical method but with different structure parameters, 
i.e. those from Table I corresponding to the columns with the exact solutions. The quoted 
parameters define a classical potential, given in the lower panel of Fig. 18, which is very 
different from the one used in the classical picture. The discrepancy is caused by the high 
anharmonicities involved. Actually, the two pictures, semi-classical and quantal, agree with 
each other only in the harmonic limit. Comparing the potentials for the two isotopes of 
Gd one finds a qualitative explanation for the behavior of the ratio i?4+ / E2+ which suggests 
that ^^^Gd is a good candidate for the critical point in the shape phase transition which 
takes place in the chain of Gd even isotopes. 

It is worth noticing that for a long time, theoretical works were focused on explaining the 
high spin states in the ground band, but not so much was done about bands aside the ground 
state band. Now we are confronted with a new situation. Indeed, to explain consistently 
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very many excitation energies of states with low angular momenta is a real challenge for 
any theoretical approach. For example in ^^^Er a large number (105) of energy levels are 
experimentally known, 26 of spin and 79 of spin 2. 

Using a sixth-order boson Hamiltonian we derived analytical formulas for the excitation 
energies of these states which involve a small number of parameters: four in the semi- 
classical treatment and five in the boson description. Both sets of formulas are describing 
quantitatively quite well the existent data. In order to draw a conclusion about how these 
coefficients depend on the atomic mass a richer systematics is necessary. 

One may argue that for many of the states considered here, the single particle degrees of 
freedom prevail. Actually we may share this opinion but, on the other hand, we think that 
the single particle behavior may be simulated by the anharmonicities involved in the present 
phenomenological model. Some of the considered states may have collective features. It is 
worth mentioning that the present model is able to account for these properties shown by 
a deformed nucleus such as ^^^Er despite the fact that one uses a boson number conserving 
Hamiltonian. Our attempt is not singular in this respect. Indeed, this is one of the signatures 
of the interacting boson approximation [19] which is successful in describing rotational bands 
in non-spherical nuclei. Moreover, the Hamiltonian given by Eq. (1) with F = 0, has been 
previously used to describe the yrast bands in transitional and deformed nuclei 20|, l21 |. 
The results of the quoted papers show that some properties determined by the nuclear 
deformation can be described by a suitable choice of the structure coefficients multiplying 
the anharmonic terms. Certainly, data concerning the electromagnetic transitions of these 
states are necessary in order to have an additional test and a more complete picture. 

A very nice test of the predictive power of our simple formulas was obtained by applying 
them to the newly found data for ^^'^Er, by keeping the numerical values for the structure 
coefficients as obtained in our previous calculations. We showed that the new data are 
surprisingly well described by the same parameters set. 

Also, very simple formulas for the B(E2) values characterizing the transitions between 
the states are derived within the two approaches. Note that the expressions for the semi- 
classical transition \2{n — 1)) \0n) and the boson transition \n — 1 102) —>■ \n 000) are 
identical. 

Of course the microscopic descriptions has the great merit of interpreting the data in 
terms of the single particle motion and may address some issues which are complementary 
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to those accessible for phenomenological models. On the other hand in phenomenologi- 
cal pictures one may find a way to improve the microscopic description. Indeed, in our 
formalism we have seen that the sixth-order boson term is necessary in order to obtain a 
quantitative description of the data. On the other hand an RPA treatment of a Hamiltonian 
involving the mean field and a two body Q.Q interaction, yields a quadratic boson term. 
Going beyond QRPA by a boson expansion procedure, higher order terms in bosons are 
obtainable. It is well known that the quadrupole two quasiparticle operator A\^{a, h) can be 
expressed as an odd powers expansion in bosons while the quadrupole quasiparticle density 
operator B2^{a,h) as an even power. Therefore, ignoring the quasiparticle correlations due 
to the operators Ag^Sl^, A2^i?2/i, ^2^1-82^, ^2^-82^1 and considering the first-order boson ex- 
pansion for the remaining terms, one obtains a sixth-order boson Hamiltonian similar to the 
phenomenological Hamiltonian used in the present paper. Concluding one may assert that 
the phenomenological Hamiltonian has actually a microscopic counterpart. On the other 
hand our calculations suggest that going beyond QRPA the microscopic model QPM might 
describe in a better quantitative way the existent experimental data. 

Recently the states 0^ have been considered by the Interacting Boson Approximation 
(IBA) approach within the phase transition context and several analytical results have been 
derived ^]. Thus it was shown that in each of the symmetries ?7(5), 0(6) and SU{?i) 
the energies of the state 0^ depend linearly on n in the regime of large values for the 
quadrupole boson numbers. Several degeneracies of the states 0^ with the states Ji have 
been pointed out for the critical values of the ordering parameters. In particular, the degen- 
eracy -E'(O^) ~ E{Qi) might be viewed as a hallmark of the X(5) symmetry and moreover 
the ratio i?(6]^)/-E(0^) could play the role of the order parameter for the specific phase tran- 
sition. For the symmetries characterizing the critical points of phase transitions, analytical 
expressions for the energies E{Q^) have been derived. The energy -E(O^) depends on n as 
n{n + x) with x depending on the considered phase transition. 



In order to make a fair comparison of the present results with those of Ref. 3a], analytical 
results of -^(0+) for any interaction strength of the IBA Hamiltonian would be desirable 
and moreover the same order boson Hamiltonian to be considered. By considering in the 
expressions 12.191 the situation when the six order term is missing, i.e. F=0, one obtains: 
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En,o,o 

En,3,0 



BN{N + x), N = 0,2,4, 

9 
-( 

4 



(5.1) 



BN{N + x) + -C,N = 3, 5, 7, .... 



with x={A + 7 - |C)/5. Thus, one may say that the results of Ref . [35|] may be recovered 
by the present formahsm in the hmit of F = 0. However, the six order term determines 
an additional minimum in the classical potential energy and therefore a new phase of the 
nuclear system is expected. Moreover, the cubic dependence may account for the complex 
structure of the energy distribution for the states 0^ and 2^. 
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VI. APPENDIX A 



Here we give the results for the matrix elements of the transition operators factors, which 
are needed for calculating the B{E2) values. The matrix elements of (3 are: 



{Fn+l l\P\FnO 
{Fn~l l\P\FnO 

{Fn+l 2\P\Fn3 

{Fn-1 2\P\Fn3 

{Fn+l 2\P\Fnl 

{Fn~l 2\(3\Fnl 
{Fn+l l\(3\Fn2 

{Fn~l l\(3\Fn2 



In + 5 



In 



ln + 6 

'n + 6 

In — 1 
2 ' 

n 

r_ 

'72+5 



(Al) 
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With the convention of Rose for reduced matrix elements, one obtains: 



(^102||^||^OOo) 



1 



(^202! I '^2 1 1^310) 



1 



(^102! 1^1 1^202) 



(^202! 1^1 1^102) 
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